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Some possible benefits of dynamical collapse for a quantum theory of cosmogenesis are discussed. 
These are a possible long wait before creation begins, creation of energy and space, and choice of a 
particular universe out of a superposition. 
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I. INTRODUCTION 

Suppose one believes that the onset of our universe 
should be described by quantum theory. Quantum evo- 
lution famously produces a state vector which is a su- 
perposition of possibilities, but our particular universe is 
just one of those possibilities. How was the choice made? 
A dynamical collapse theory can help with this. 

Such a theory as CSL (Continuous Spontaneous Local- 
ization) constructed to describe wave function col- 
lapse within the present universe, works quite well. The 
state vector evolves continuously, governed by a Hamil- 
ton whose Hermitian part is the usual Hamiltonian, and 
whose anti-hermitian part depends upon a randomly fluc- 
tuating field. When the hermitian Hamiltonian evolves 
the state vector into a superposition of macroscopically 
different states, the anti-hermitian Hamiltonian causes a 
rapid evolution to one of them. 

A dynamical collapse theory posits that the state vec- 
tor corresponds to the real state of the world, and pro- 
vides a mechanism to accomplish this. If one believes 
that nature uses the collapse mechanism to choose the 
outcome of an everyday event, it is natural to suppose 
that a similar mechanism was at work for the Big Event. 

Suppose the pre- universe is modeled by a state |0) 
which represents no space. Further, suppose that time 
exists. (An eventual hope might be to have time come 
out of such a theory, but crawling comes before walking.) 

A Hamiltonian must be chosen to evolve the state |0) 
into our early universe. One might suppose this is a rel- 
atively rare event, i.e., the universe has a low probabil- 
ity/sec of being created out of the pre- universe. Dynam- 
ical collapse can help with this. 

In standard quantum theory with the instantaneous 
collapse postulate, there is a phenomenon called the 
"Zeno" or "watched pot" effect. Suppose a system starts 
in state |0) and evolves under Hamiltonian H for a short 
time At, followed by instantaneous collapse. Suppose 
this occurs repeatedly over the time interval T. It is 
straightforward to show that the probability the system 
remains in state |0) is exp — [(AiJ) 2 TAt + o(At) 2 ], where 
{AH) 2 ee (0\H 2 \0) - (0|ff|0) 2 . Thus, as At gets smaller, 
the probability that |0) evolves is reduced. A compara- 
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ble behavior occurs for continuous collapse. When the 
collapse rate parameter A increases sufficiently, the prob- 
ability that the initial state vector evolves is reduced. 

Not only space needs to be created, but also the energy 
which fills that space. Dynamical collapse can help with 
this. 

In usual quantum theory, energy is conserved. How- 
ever, that need not be the case in a dynamical collapse 
theory. For example, in CSL, particle wave functions are 
narrowed by the collapse process, which means that their 
energy grows [1]. 

In this paper, a simple cosmogenesis model to illustrate 
these features is presented. Its mathematics is that of 
a spin 1/2 caused to rotate with angular frequency u> 
by a constant y-directed magnetic field, and that of a 
perturbed harmonic oscillator, both subject to collapse 
dynamics. 

The spin 1/2 acts like a gatekeeper to block (state | ^)) 
or allow (state | \.)) the harmonic oscillator to operate. 
The spin starts in the block state but is prevented from 
rotating with frequency uj because of "watched pot" be- 
havior: a sufficiently large value of the collapse rate pa- 
rameter A can make the probability /sec of a transition 
to the allow state arbitrarily small. 

The harmonic oscillator starts in the ground state |0) 
representing no space volume. The states |n) represent 
a universe of n cells, with energy en and volume vn. 
To be concrete, we shall suppose that e=Planck energy. 
v=Planck volume is a possible choice, or v may be an 
increasing function of time. One may consider that the 
energy generated is the dark energy alone, or all the en- 
ergy: the nature of the energy and the size of v do not 
play a role in the dynamics (but see section IVI]) . 

Due to a perturbation, the oscillator would just os- 
cillate between |0) and slightly higher \n) states. (Two 
kinds of perturbation shall be considered, an added term 
proportional to position in section HTT1 and to squared po- 
sition in II VI ) However, because of the attendant collapse 
(rate parameter A'), the universe is generated: the state 
vector evolves states of increasing values of n. With the 
initial state | f> |0), the (unnormalized) state vector at 
time t is taken to be 

. e ~So dt'{i[£ata+g(at+a) s ]Q+ 3 i 7 [u.'(t')-2A'a t a] 2 }| ^ > /j\ 

In Eq.|T|), T is the time-ordering operator. If there 
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were no collapse, the spin would rotate clockwise in the 
x — z plane. Q = | J,) (J, | is the projection operator 
on the allow state, g is a coupling constant: we shall 
consider the cases of s = 1 (an interaction which displaces 
the harmonic oscillator) and s — 2 (an interaction which 
changes its spring constant). 

w(t') and w'(t') are white noise random functions of 
time. The probability density for w, w' to take on par- 
ticular values over the time interval (0 < t' < t) is given 
by the probability rule: 

P{w(t), w'{t)) ~ w , M /<*, t\V, t) w , w , (2) 

where the proportionality constant is determined by 

/oo 
DwDw'P = 1 

with DwDw' = Yii dw(ti)dw' (ti) (U+i — ti = dt). 

According to Eq.([T]), each different w(t'), w'(t') gives 
rise to a different state vector \^,t) w ^ w i. The probability 
of that state vector (of that w(t'), w'(t')) is given by 
Eq.@. Eqs. (|ll2[) completely define the model. It only 
remains to draw the consequences. 

In this paper, we shall only consider the evolutions 
of the spin and harmonic oscillator separately, not their 
joint evolution, in order to keep the discussion simple and 
to elucidate the features of each. 



A. uj = 

We begin by considering uncontaminated collapse be- 
havior (oj = 0). Although Eq.Q is exactly soluble (see 
5]. restated in Appendix [Bl in the notation of this pa- 
per), we prefer to examine its properties in a way that is 
transferrable to the w^O case. 

Here is the general context @. Consider an ensemble 
of normalized state vectors \ip,t) = Ylk=x c k(t)\k) which 
evolve stochastically, but start out identically at t = 0. 
Two conditions on Xk(t) = \ck(t)\ 2 are sufficient to obtain 
collapse behavior obe ying t he Born rule. These are the 
Martingale condition Xk{t) = Xk{0) and the asymptotic 
correlation condition a;fc(oo)a;;(oo) = [k ^l). 

The latter equation is an integral over the product of 
non-negative terms. The only way it can be satisfied is if, 
for every state, every Xkipo) but one vanishes. Therefore, 
the asymptotic probability density distribution must be 

P{xi, ...x n ,oo) = aiS(l - xi)S(x 2 )...S(x N ) + ... 

The Martingale condition then implies — Xfc(0), i.e., 
the Born rule. 

For our problem, there are two states: x\(t) = 
cos 2 [#(i)/2] is the squared amplitude of | |) and x 2 {t) — 
sin 2 [6>(£)/2] is the squared amplitude of | |). We use 

(f4]) to calculate cos8(t), sin 2 9(t). After integrations by 
parts, we obtain: 



II. GATEKEEPER/SPIN 

A numerical analysis has been given for this 
problem 4]. Here we give an analytical treatment. 

The Fokker-Planck equation for the probability distri- 
bution R(9,t) of the angle < 9 < 2ir that the spin 
makes with the z-axis is 



d_ 

d9 



f«(M) 



\gp_ 

2W 



(dey 

dt 



R(6,t) 



(3) 

where the average is over all possible trajectories of 9 due 
to all possible w(t). Appendix |A1 shows how to calculate 
the drift d6/dt and the diffusion dfP/dt (Eqs.((ME3), 
yielding the Foker-Planck equation: 

^-R = [u - 2A sin 9 cos 6]R + 2X-^-r sin 2 9R. (4) 
dt oO o9 z 

Qualitatively, we see that the drift's steady clock- 
wise rotation with angular velocity uj is opposed by the 
stochastic drift in the first and third quadrants, and aided 
in the second and fourth quadrants: the stochastic drift 
tends to drive the spin to the nearest pole (9 — 0,tt). 
The diffusion rate is largest at z = 0. The diffusion rate 
approaches at the poles ~ 9 2 , (tt — 9) 2 . This makes the 
poles absorbing boundaries [(|: if uj = 0, all the trajecto- 
ries eventually end up at the poles (collapse). 



d d f 2 * 

—cos9(t) = — d9 cos 9R = 
dt w dt J 

8 8 f 2 * 

— sin 2 9(t) = — / d9 sin 2 9R 
dt w dt J 



(5) 



-4Asin 4 (9(t). (6) 



It follows from Eq.([5]) that the squared amplitude 
(1/2)(1 + cos (9) = cos 2 (#/2) is constant, so the Martin- 
gale condition is satisfied. 

Now, sin 2 8{t) is non- negative but, according to Eq.(j6]), 

it keeps decreasing unless sin 4 9(t) vanishes. Therefore, 

sin 4 9(t) must vanish at some time t = T. But, that is 
only possible if R(9, T) is non-zero only where sin 9 is 

zero, i.e., at 9 = 0,7r. But, this means that sin 2 9{T) = 
too. Since the average of the product of the up and down 
amplitudes is sin 2 [6>/2] cos 2 [(9/2](T) = (l/4)sin 2 9{T) = 
0, the asymptotic correlation condition is therefore satis- 
fied. Thus, there is collapse satisfying the Born condition. 



B. uj / 



It follows from Eq. (JH) that 
d 



dt 

m 



cos 6{t) 



ujsn\8{t) 



sm9(t) = ujcos9(t) - 2Asin0(t), 



(7) 

(8) 
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whose solution is 



cos8(t) = - — 3_ — | cos6>o[A+e 
— A_ I 



A_e 



-x+u 



-ui sin #0 e 



A± = \ ± \/\ 2 - uj 2 



-A +t] | 5 



(9) 



2A, — for A >> uj. 
2A 



where 0q is the value of 8 at t = 0. 

Consider the initial state 9q — 0, spin up. The ensem- 
ble average squared amplitude of | f) is cos 2 [0/2] (i) = 
.5[1 + cos 6(t)]. Thus, according to Eq.©, 



cos 2 [6>/2] (t) = 



2(A+-A_) 



{a+[i 



,-A-tl 



A_[l 



,-Aj.ti 



1 + e" 



-(w 2 /2A)t 



(10) 



Eq. (fT0| shows the "watched pot" behavior mentioned 
earlier. For ui 2 /\ « 1, most trajectories remain in the 
neighborhood of = 0, i.e., there is only a small proba- 
bility for a transition from | f ) to | T) < 

Although it is an irrelevant feature for our application 
of this model, we close by remarking that, eventually, 
an equilibrium distribution of spins is reached, which is 
independent of the original state. When R(9, t) is taken 
to be independent of i, Eq.Q provides the first order 
differential equation 

Q 

- [u - 2A sin 9 cos 9]R eq (6) + 2A— sin 2 0R eq (6) = -C, 

do 

(11) 

</:- — — r (12) 



whose solution is 



e 2\' 



[l + (z-cot0) 2 ] 3/2 



According to Eq. (ll2[) . the equilibrium probability distri- 
bution is largest at 6 — 0+, tt+, and falls off as 6 increases 
from to it— and from 7r+ to 2ir— , with a step- up dis- 
continuity crossing the poles in the clockwise direction. 



III. SPACEGENERATOR/OSCILLATOR I 

The state vector evolution we consider here is, from 
Eq.©, 

= re" /o* ««'{*»+ AK*0-2AJV] a }|o). (13) 

where H = eN + g(a^ + a) and N = a^a. The density 
matrix follows from Eqs. (lll2l) : 



p{t) = / Dw w (*,t\*,t), 



\^>,t) ww (^,t\ 

f e - fo dt'{i{H L -H R ] + ±({N L -N R ] 2 }\ Q ^ Q ^ ( 14 ) 



where an operator with subscript L (R) acts to the left 
(right) of p, and T time reverse-orders operators to the 
right of p. Therefore, from (|14[) . the density matrix evo- 
lution equation is 

^!± = -i[H,p(t)}-~[N,[N,p(t)]\. (15) 

With the notation (A)(t) = Tr&ce[Ap(t)], we find the 
coupled equations 



d(N)(t) 



dt 



d(a)(t) 



ig[(a)(t) - <at)(t)] 
A- 



(16) 



(U = -2e(a)(t)- ig --(a)(t). (17) 

Solving Eqs.dTHl fTTj) and the complex conjugate of (jTT|) 
results in: 

2ff 2 fA 



{N){t) = (A/2) 2 + ( 2 {f- cosB + e ^ l cos[et-e]} (18) 

where tan# = e/A. 

We see from (fT8"|) that, when A = 0, there is only os- 
cillation, but when A / there is linear growth of the 
ensemble mean of the expectation value of the energy eN 
and the volume vN. The mechanism is that the coupling 
~ g causes N to oscillate by a small amount away from 
its initial value 0, but there is a chance that collapse will 
occur to an N ^ value, and then the oscillation takes 
it away from that, etc. (It is worth noting here that 
"watched pot" behavior is again in evidence, i.e., A — > oo 
implies JW) -> 0.) 



IV. SPACEGENERATOR/OSCILLATOR II 

In the model in the previous section, the mean vol- 
ume of space grows linearly with time. Although infla- 
tion models make space expand exponentially with time, 
they do not start from zero space as do the models un- 
der consideration here. Inflation models assume that 
enough space had been created earlier so that it may 
subsequently be described by classical general relativity, 
with the inflaton field which lives on it obeying quantum 
field theory. 

There is nothing wrong with an initial space creation 
that grows linearly with time. Nonetheless, it is inter- 
esting to see that an exponential growth can be obtained 
with only a slight modification of the model of section \TU\ 
replace the perturbation in the Hamiltonian by g(a^+a) 2 . 

One then obtains from Eq. ([T4f the coupled equations 

d W) o...r7I3V 



dt 

dja 2 ) 
dt 



2ig[(a 2 ) - (at 2 )] (19) 
2ja 2 )[ie + 2ig + 2A] - 2ig[l + 2(iV)] (20) 
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and the complex conjugate of (|20|) . Writing all three 
quantities as ~ exp st and putting that into the homoge- 
neous part of the equations results in the cubic equation: 

/(s) ee s 3 + 4As 2 + (e 2 + A 2 + Aeg)s - 32A 5 2 = 0. (21) 

For example, if g = 3A/2 and e = [y/U - 3] A, Eq.flHl) 
becomes 

[s - 2A][s - 6Xe l27v/3 ][(s - 6Ae~ l2T/3 ] = 0. 

For another example, if A << g = e/4, then there is a 
root for which the s 3 and s 2 terms may be neglected, with 
the remaining terms giving s sa 32A<? 2 /(e 2 + 4eg) = A. 

Generally, it follows from Eq. (j2"Tl) that there are two 
complex conjugate roots with negative real parts. Be- 
cause /(0) = — 32A<7 2 < and / has positive slope for 
s > 0, there is one solution of (|2"Tj) which is a positive real 
root, giving rise to an exponential growth of the mean en- 
ergy and volume of space. 

V. ON CHOOSING THE UNIVERSE 

As an indication of how the collapse acts to " choose a 
universe," in either model, suppose we neglect the action 
of the perturbation ~ g over the interval (t,to). Then, 
the off-diagonal components of the density matrix decay 
according to the solution of Eq. (TT5|) , 

(n\p{t)\m) = e~ (t - to)[le(n ^ m)+ *(™- m ) 2 ](n|p(t )|m). 

(22) 

Of course, the density matrix approaching diagonal form 
does not guarantee that each individual state vector in 
the ensemble is approaching some definite state |n), but 
that is the case hereQ- 

However, one cannot neglect the action of the pertur- 
bation: it is working alongside of the collapse. As we 
have seen in section [21 when collapse competes with the 
Hamiltonian evolution, the outcome is not obvious. 

What is needed if the collapse is to "choose a uni- 
verse?" . A universe is a macroscopic object and, as such, 
should have well-defined values for macroscopic variables 
such as the amount of mass it contains. That is, for any 
probable state vector in the ensemble, although it is a 
superposition of states of various n-values, its range of 
n-values should be small compared to its mean n-value. 
Typically, this requirement is quantified by requiring the 
deviation from the mean to be small compared to the 
mean. 

Let's state this requirement in the context of our prob- 
lem. For any operator A, define its expectation value as 

(A)(t) = w (y,t\A\V,t) w / w (<l>,t\y,t) w . 

Define the ensemble average of products of such expec- 
tation values as 



(A)(B)... ee / Dw w (V,t\*,t} w (A}(B).... 



Then, the squared standard deviation of A is 

(AN) 2 ee ((A - (AO) 2 ) = W 2 ) ~ 

and the requirement is AN/(N) << 1. 

It should be emphasized that (AA) 2 is not 

(A A) 2 ee ((A -(AO) 2 ) = W) ~ W) 2 - 
Indeed, (AN) 2 < (A A) 2 , since 

< [(A -"(A)")] 2 = (A) 2 - "(AO 2 . 

As we have seen, using the density matrix it is easy to 
calculate ensemble averages involving two state vectors. 
However, it is not easy to calculate ensemble averages 
when four or more state vectors are involved. 

In the two previous sections we have sho wn h ow to 
calculate (N)(t), and can similarly calculate (N 2 )(t), so 
(AoA) 2 is easy to find. It turns out that, for large t, for 

both models, (A A) 2 ~ (JV) (i.e., ~ t 2 for the first and 
~ exp 2st for the second) . 

However, (N) 2 (t) involves the ensemble average of a 
product of four state vectors, so (AA) 2 is not easy to 
find. 

Here is what we have been able to do. In Appendix [Cl 
it is shown in Eq. (|C9|) that 

i^ 2 - >^ _ 

+2» ([A-(A),g])(A-(A)) 
-4A((A-(A0) 2 )((A-(A0) 2 ). (23) 

This would seem to increase the difficulty since, e.g., the 
last term in (|23p is the sum of terms involving the average 
of products of four, six and eight state vectors. However 
let us make the not unreasonable assumption that, for 
most state vectors of high probability, (A — (A)) is nearly 
constant, and likewise for ((A — (A)) 2 ). Then 

(A)(N-jN))^jA)(N-jN))=0, 

((A-(A)) 2 )((A-(A))2) « [((A — (A)) 2 )] 2 , 
and Eq. (j2"3")) becomes 

^(AA) 2 = ^-(A A) 2 - 4A[(AA) 2 ] 2 . (24) 
at at 

This is an example of the Riccati equation dz / dt + z 2 = 
fit). It may be converted to a linear second order equa- 
tion by the transformation z = ii/u, its solution for our 
two examples for large t involves Bessel functions, etc. 
However, the important point is that, for our two ex- 
amples for large t, the left side of Eq.(|24| is negligible 
compared to either of the two terms on the right side, so 

^Sv^ ^ (25) 
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For the example in section IIII1 we may write, 
{A Q N) 2 pa (N t/T ) 2 for large t, where N Q is the cur- 
rent number of Planck masses in the universe and To is 
the current age of the universe. We therefore find from 
Ed.([25)l that 

(N) [2XT Q ]^N 1/2 

For the example in section IIV1 we may write, 
(A A^) 2 fa [N a exps (t - ?o)] 2 for large t, where s N 
is the current rate of increase of Planck masses (see the 
next section for more discussion of this). We therefore 
find from Eg. ([25)1 that 



VW^ (N) [2A/s ]W 1/2 ' 

As mentioned in the next section, it is plausible 
that XT and X/sq are of order 1. Therefore, indeed, 

AN/jN)^ No 1/2 « 1. 

VI. CONCLUDING REMARKS 

By means of a simple quantum mechanical model, the 
usefulness of collapse dynamics in creating the universe 
out of nothing after a long wait, generating its energy, 
and choosing a macroscopic universe out of a superpo- 
sition, has been sketched. The model may be regarded 
as just an illustration, but one might try to take it more 
seriously. 

It is perhaps natural to choose A -1 and g" 1 to be of 
the order of the Planck time T p . If one takes the volume 
of the present universe to be Vq ~ 3 x 10 80 cubic meters, 
takes the mass-energy density to be the present critical 
density value po pa 10~ 26 kg/m 3 then, with the Planck 
mass e = M p pa 2 x 10~ 8 kg, there are presently No = 
poVo/Mp pa 10 62 Planck masses in the universe. 

Suppose that all this mass-energy is generated at the 
beginning of the universe, according to the model of sec- 
tion HVl N(t) pa exp(*/T p ). It would take time t sa 140T p 
to generate Nq such masses. If each Planck mass occu- 
pies a Planck volume, the universe at that time would 
be of radius pa (3iV /47r) 1 /3£ p ~ 5 x l(r 13 cm (L p Pa 
1.6 x 10 -33 cm is the Planck length), about 5 times larger 
than the radius of a proton. Of course, with this scenario, 
one wants the process to abruptly stop at this time. One 
might make g, X and e proportional to Q(Nq — N), and 
let a different Hamiltonian take over the subsequent evo- 
lution but this seems rather artificial. 

More interesting would be if the process were ongoing, 
with A a suitably decreasing function of N, to provide a 
mechanism which continually slows the energy produc- 
tion down, yielding the correct amount of energy and 
energy production rate at the present era. 

With regard to the current energy production rate, 
here is a suggestive result. 



According to the model, given the present time To, 
the number of Planck masses at time t > To is N(t) = 
A r oexpAo(t — Tq), where Ao is the current collapse rate. 
Thus, since the energy in the universe is E(t) = M p N(t), 
its present rate of change is 

dE/dt = E X . 

On the other hand, suppose next that the radius of the 
universe is expanding with the speed of light, dR/dt = 
c. Then, the present rate of increase of energy in the 
universe is 

dE/dt = p dV/dt = poAnRlc = 3E c/R . 

Equating these two expressions for dE/dt, using Rq pa 
(3V0/471-) 1 / 3 = 4 x 10 26 m, we obtain 

A = 3c/i? « 2 x 10~ 18 sec"\ 

approximately the inverse of the age of the universe. 

This is not far from the value A sa lCP 16 sec _1 sug- 
gested for the GRW theory Q and adopted for CSL. It 
suggests considering the consequences if the collapse rate 
in CSL, hitherto considered a constant is, instead, the 
inverse of the age of the universe. This would imply a 
large rate of excitation of matter at the beginning of the 
universe. However, it would not have much effect, either 
on visible excitation or on the collapse of macroscopic 
objects over a time interval, say, 1/100 to 100 times the 
present age of the universe. 

The reader will have noted that there has been no men- 
tion of general relativity. Connection with the FRW for- 
malism can be made through identifying an operator for 
the "radius R" of the universe. 

R depends not only on the operator N(t) but also on 
the operator representing the volume v(t) of a cell, v 
must have dynamics because, if it equals the Planck vol- 
ume at the beginning of the universe, it does not re- 
main so. The current value of v may be found as fol- 
lows. The change of V is dV = vdN, while the change 
of energy is dE = edN. Since dE/dV = po, we obtain 
11 = e/po « 2 x 10 18 m 3 , i.e., cells of Planck energy which 
are currently added have length sa 1000km. 

So, there must be an addition to the Hamiltonian to 
give dynamics for the operator v, presumably guided by 
the FRW formalism, with 

4tt o, , /•* ,dN(t') , , s 

Moreover, if the created energy is to include not only dark 
energy but also visible matter and dark matter, dynamics 
must be included describing how the latter two evolve 
from the total energy Ne. 
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Appendix A: Spin Problem 



Appendix B: Solution of Eq.(|4]) for cu — 



Consider an initially normalized state vector \^f,t) de- 
scribing a spin in the x — z plane. If < 9 < 2ir is the 
angle between the z— axis and the spin vector, the compo- 
nents of the spin in the a z basis are (f \&,t) — cos(#/2), 
(4- t) = sin(t9/2) (up to a phase factor). The spin part 
of the evolution given in Eq. ([T]), for an infinitesimal time 
interval dt, is: 

\%,t + dt) w = e -dt{^2+M^)-^?}\^ >t ) w . (Al) 



Eq. (|Al|) . expressed in the <j z basis, retaining all terms 
necessary for first order in dt is: 

Acos[0(t + dt)/2] = e-^ w2(t) {cos(6/2)[l + w(t)dt 
+-(w(t)dt) 2 - Xdt] - sm(6/2)[-u)dt]} 

Asin[6(t + dt)/2] = e-^ w2w {sm(9/2)[l - w{t)dt 

+ ^(w(t)dt) 2 - Xdt] + oos(d/2)[io;<ft]}(A2) 

The amplitude A is necessary because the state vector 
evolution is not unitary. Indeed, this is what makes the 
probability rule Eq. ([5} meaningful: 

P(w(t))~ w {$,t\V,t) w = A 2 

= e - ix«' 2 ( t ){i + 2w(t)dt[cos 2 (9/2) - sin 2 (6»/2)] 

+2[-(w(t)dt) 2 - Xdt] + (w(t)dt) 2 } 

= e -^ w2{t) [l + 2w{t)dt cos 0] (A3) 

We note that the total probability is 1, as may be seen 
by multiplying Eq. ([AH by (2ttX/ dt)~ 1/2 dw(t) and inte- 
grating over (—00,00). In the last step of Eq. (| A3|) . we 
have set (w(t)dt) 2 = Xdt, since that is what its integral 
gives to order dt in all of our calculations (the same thing 
could have been done in Eq. (| A2|) and shall be done in 
subsequent equations). 

Solving Eqs. (|A2I) to order dt yields: 

6{t + dt) = 2tan~ 1 {tan(6>/2) + ^ujdt sec 2 (9/2) 

+ 2[Xdt-w(t)dt]tan(9/2)} (A4) 

Using tan _1 (tan <fi + x) — <fi + xcos 2 <fi—x 2 sin^cos 3 <f>+... 
in Eq. (|A4I) results in: 

d6(t) = udt - 2w(t)dt sin 9 + 2Xdt sin 9 cos 9 (A5) 



Now that we have Eqs. (|A3l[A5|) . we can find d6, (d9) 2 , 
which are used in constructing the Fokker-Planck equa- 
tion for the probability distribution of 9 in section HI1 



d9 



1 



(dey 



00 y/2nX/dt 

DO J 



— OO y/2TlX/dt 



dw(t)P(w(t))d9(t) 

= ujdt - 2Xdt sin 9 cos 9 (A6) 

dwP(w(t))(d9) 2 (t) 



= AXdt sin 2 



(A7) 



For convenience, we give here the solution^ of 

^R = 2x{^[sin9cos9}R+-^ sin 2 9R~j. (Bl) 

By changing variables to the squared amplitude x — 
cos 2 (9/2), and R = (1/2) sin9R', the equivalent equation 
is obtained: 



!# = 8A|^(l-:r)] 2 #. 



(B2) 



where < x < 1. Eq. (|B2l) clearly shows the no-drift 
condition dx(t)/dt — necessary for Born-rule-obeying 
collapse, as discussed in section Hi A[ It also shows how 
the diffusion vanishes at x = 0, 1, providing absorbing 
boundary conditions. 

A further change of variables, z = ln[a;/(l — x)] and 
R' = x(l — x)R" provides an equation with normal dif- 
fusion, 

Ft R " = 8X {& R " ~ ITz tanh (^/2)i?"}, (B3) 
which is readily solved: 



R"(z) = 



1 



z o/2 e 35Xt[ Z -z -8At] 2 



8V27rAicosh(z /2) 



We see from Eq.([B4j that R"(0) S(z - z ) and 

R"(oo) -)■ r^- , 

v ; cosh(z /2) 

In terms of 9, the solution (fB4l) of Eg .([EI]) is 

1 



J e z °/ 2 5(z - 00) + e- z °/ 2 S(z + 00)} . 



R(9)d9 

|cos 2 (6l /2)e« 



=dlntan(6»/2) 



8ttA 

— 1 r 1 t 

2//3_ SXt [ tan(0 /2 



+4At 



+ sin 2 (6lo/2)e^ 







.if),,/ 2 



-4At 



'}.(B5) 



We see from Eq. (jB5]) that R(0) -> 8(9 - 9 ) and note the 
collapse behavior obeying the Born rule, 

R"(oo) -> cos 2 (6»o/2)(5(6i) + sin 2 (6' /2)5(6» - tt). 



Appendix C: Ensemble Averages of Products 

We first prove the following theorem. Consider the 
state vector evolution (TT3l). Then 



dt 



(A)(B) = (-i[A,H] (X/2)[N, [N,A]])(B) 
+ (A)(-t[B,H]-(X/2)[N, [N, B}]) 



-AX (N) 2 (A)(B) 



-2X[(N) (A) (BN + NB) + (N) (B) (AN + NA)} 
-X[(AN + NA) (BN + NB)}. (CI) 
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InEq.dUT]), (Q)(t) = w (^,t\Q\^!, t) w / w (% t\*,t) w . The 
overline refers to the ensemble average: for example, 
(3)(t) = f Dw w ($,t\V,t) w {A)(t) =1tAp(t), where p(t) 
is the density matrix. 

According to Eq. (fl3l) , the state vector at time t + dt is 



Writing the last bracket in Eq. (|C6|) as 



+4A*(A^) 2 (t)(C7) 



e 4A 



■[1 -idtH + dtw'N- (X/2)dtN 2 ]\^,t) w . (C2) 

In Eq. (|C2|) , we have written zu' as the scalar field at time 
t + dt, and w represents the scalar field for < t < t. We 
have also set (w'dtN) 2 /2 = XdtN 2 /2, and will continue 
to do so, since that is what it will equal when integrated 
over the probability. 

The ensemble average of the expectation value of .4 
at time t + dt (calculated using Eq. (|C2[) ). is obtained by 
multiplying the expectation value, 

(A)(t + dt) = (A)(t) + dt{-i([A,H])(t)+w'({N,A}+){t) 
-(X/2)([N, [N,A(t)]])(t)} 

r w {V,t\$,t) w -I 



the integrals can be performed with the result Eq. flCip . 

One may check that, if B = 1 is inserted into Eq. (|Cll) . 
then Eq. (|C5j) is obtained. 



Now, set A = B = N. Eq.lKJil becomes: 



dt 



(N)i(t) = -2i([N,H])(N) 



+4A((A — (N)) 2 )((N — (A)) 2 



I to {*,t|(l + 2dhi;'J\0|*,i)t 
by the probability density 



(C3) 



Since 



(y,t + dt\y,t + dt} w i iW 

2dtw' N -2\N 



= P. 2X 



V,t\e zatw "~ ZA "~\^,t) w 
V,t\(l + 2dtw'N)\V,t) w (C4) 



((N - (N)) 2 ) = (N 2 ) - (N) 2 



and by dw' Dw, and integrating over all w' , w, to obtain 
the well-known result: 

d- 



and 



dt (A)(t) = Tr[{-i[A,H] - (A/2) [A, [N,A(t)]]}p(t)p5) 

Similarly, the ensemble average of the product of two 
density matrices at time t + dt is 



((N - (N))*) = <iV 2 ) - (N) , 



it follows from Eq.(|C8|) that 



(A)(t + dt)(B)(t + dt) = J dw'Dwe w 2X (*,t|*,t) 
■{(A)(t) + dt{-i([A,H])(t) + w'({N,A} + )(t) 

-(\/2){[N, [N,A(t)]})(t)}} 
■{(B)(t) + dt{-i([B,H])(t) + w'({N, B}+)(t) 

-(\/2){[N, [N,B(t)]])(t)}} 

l w {$,t\(l-2dtw>N)\V,t) w y 



-((TV- (AO) 2 ) = -((TV- (AO) 2 ) 



(C6) which is used in section V, Eq. (|23|) . 



(C8) 



-2i ([N-(N),H])((N-(N)) 
-4A<(JV - <JV)) a X(JV - <JV» a >- (C9) 
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